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Abstract 

By analyzing the survival probability amplitude of an unstable 
state we show that the energy corrections to this state in the long 
(t — ► oo) and relatively short (lifetime of the state) time regions, are 
different. It is shown that in the considered model the above correc- 
tions decrease to zero as t — > oo. It is hypothesized that this property 
could be detected by analyzing the spectra of distant astrophysical 
objects. The above property of unstable states may influence the 
measured values of possible deviations of the fine structure constant 
a as well as other astrophysical and cosmological parameters. 

PACS numbers: 03.65.-w, 11. 10. St, 06.20. Jr, 98.62.Py 
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1 Introduction 

Within the quantum theory the state vector at time t, \<f>(t)), for the physical 
system under consideration which initially (at t = to = 0) was in the state 
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|$) can be found by solving the Schodinger equation 

ift^|$(t)> = J ff|$(t)) > |$(0)) = |$), (1) 

where |$(t)), |$) € TC, H is the Hilbert space of states of the considered sys- 
tem, || |<&(t)) || = || |$)|| = 1 and H denotes the total selfadjoint Hamiltonian 
for the system. If one considers an unstable state |<3>) = \4>) of the system 
then using the solution | </>(£)) of Eq. (JTj) for the initial condition 1 0(0) ) = \<p) 
one can determine the decay law, V</>(t) of this state decaying in vacuum 

V 4l {t) = \a{t)\\ (2) 

where a(t) is probability amplitude of finding the system at the time t in the 
initial state \<p) prepared at time to = 0, 

a(t) = (ct>\<f>(t)). (3) 

We have 

a(0) = 1. (4) 

From basic principles of quantum theory it is known that the amplitude 
a(t), and thus the decay law V^it) of the unstable state are completely 
determined by the density of the energy distribution u(S) for the system in 
this state pQ, 

a{t) = [ u(£) e~i £t d£. (5) 

JSpec.{H) 

where u>(£) > 0. 

Note that (jSJ) and (jlj) mean that there must be 



o(0) = / u(£)d£ = 1. (6) 

From the last property and from the Riemann-Lebesgue Lemma it follows 
that the amplitude a(t), being the Fourier transform of tu(£) (see (jSD), must 
tend to zero as t — > oo pQ. 

In [2] assuming that the spectrum of H must be bounded from below, 
(Spec.(H) > — oo), and using the Paley-Wiener Theorem [3] it was proved 
that in the case of unstable states there must be 

\a(t)\ > Ae~ ht \ (7) 
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for \t\ — > oo. Here A > 0, 6 > and < q < 1. This means that the decay 
law "P</,(t) of unstable states decaying in the vacuum, d2J), can not be described 
by an exponential function of time t if time t is suitably long, t — > oo, and 
that for these lengths of time V<j,(t) tends to zero as t — > oo more slowly 
than any exponential function of t. The analysis of the models of the decay 

processes shows that V^{t) — e ft , (where 7° is the decay rate of the state 
|0)), to an very high accuracy for a wide time range t: From t suitably latter 

7° 

than some T ~ to = but T > to up to t — -4- and smaller than 

t = T^, where denotes the time t for which the nonexponential deviations 
of a(t) begin to dominate (see eg., [2], [3] [8]). From this analysis it follows 
that in the general case the decay law V^t) takes the inverse power-like form 
t~ x , (where A > 0), for suitably large t > 3> t$ [2], [1] [7]. This effect is 
in agreement with the general result ([7]). Effects of this type are sometimes 
called the "Khalfin effect" (see eg. 0). 

Note that one can force V${t) to have an exponential form for all t in- 
cluding t — > 00 by using the time dependent decay rate [91 [10] . 

Unfortunately, there is no experimental evidence for the Khalfin effect. 
The performed test of the decay laws did not indicate deviations from the 
exponential form of V^>(t) at the long time region [HI [121 IB], neither did they 
show a time dependence of the decay rate at this time region. Following 
the estimations given in [5] one infers that this probably occurs because it 
is almost impossible to build in a laboratory a device detecting decays of 
the unstable states in the vacuum which allows one to observe these decays 
from the instant t = t = of the preparation of these states up to times 
t ~ 10 n Trf,, where n > 9. 

Note that in fact the amplitude a(t) contains information about the decay 
law V(f,(t) of the state \<p), that is about the decay rate 7° of this state, as 
well as the energy £® of the system in this state. This information can be 
extracted from a(t). Indeed if \<p) is an unstable (a quasi-stationary) state 
then 

a(t) ^ e~^ 8 l ~ 2^) * (8) 



So, there is 



2 dt a(t) 

in the case of quasi-stationary states. 



po_l^o = in 5a(t) J_ 
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The standard interpretation and understanding of the quantum theory 
and the related construction of our measuring devices are such that detect- 
ing the energy £® and decay rate 7° one is sure that the amplitude a(t) has 
the form (jSj) and thus that the relation (Q occurs. Taking the above into ac- 
count one can define the "effective Hamiltonian" , h^, for the one-dimensional 
subspace of states H\\ spanned by the normalized vector |0) as follows (see, 
eg- IW\) 

, def . fc da(t) 1 

h * = th -dTW (10) 

In general, can depend on time t, = h$(t). One meets this effective 
Hamiltonian when one starts with the Schrodinger Equation (pQ) for the total 
state space 7i and looks for the rigorous evolution equation for the distin- 
guished subspace of states Tt\\ cW. In the case of one-dimensional TC\\ this 
rigorous Schrodinger-like evolution equation has the following form for the 
initial condition a(0) = 1 (see [10J and references one finds therein), 

jftMfi = h+{t)a{t). (11) 

Relations (flOl) and (ffTT) establish a direct connection between the amplitude 
a(t) for the state |</>) and the exact effective Hamiltonian h<j,(t) governing the 
time evolution in the one-dimensional subspace Ti.u 3 Thus the use of 
the evolution equation fflTj) or the relation fflOl) is one of the most effective 
tools for the accurate analysis of the early- as well as the long-time properties 
of the energy and decay rate of a given qausistationary state \4>(t)). 

So let us assume that we know the amplitude a(t). Then starting with 
this a(t) and using the expression ffTUl) one can calculate the effective Hamil- 
tonian h$(t) in a general case for every t. Thus, one finds the following 
expressions for the energy and the decay rate of the system in the state \4>) 
under considerations, 

= ^(t) = »(^(t), (12) 
70 = 7*(*) = -23(M*)> (13) 

where 3? (z) and S (z) denote the real and imaginary parts of z respectively. 

As it was mentioned above the deviations of the decay law V^{t) from the 
exponential form can be described equivalently using time-dependent decay 
rate. In terms of such 7^(t) the Khalfin observation that V^t) must tend 
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to zero as t — > oo more slowly than any exponential function means that 
70 (t) < 7° for t > ^ and lim^^ 7^) = 0. 

The aim of this note is to examine the long time behaviour of £<p(t) 
and y$(t) using a(t) calculated for the given density u){8). We show that 
£tj,(t) — > as t — > oo for the model considered and that a wide class of 
models has similar long time properties: ^(t)!^^ ^ 8®. It seems that in 
contrast to the standard Khalfin effect [2] in the case of the quasistationary 
states belonging to the same class as excited atomic levels, this long time 
properties of the energy £$(t) have a chance to be detected by analyzing 
spectra of very distant stars. 

The paper is organized as follows. In Sec. 2 the above mentioned model 
is considered. Sec. 3 contains an analysis of some general case of a(t) causing 
£(j,(t) to tend to zero as t — > oo. A discussion and final remarks can be found 
in Sec. 4. 



2 The model 

Let us assume that Spec.(H) = [0, oo) and let us choose o>(£) as follows 



o 



N 1* 



u(S) = - 6(8) " , (14) 

where is a normalization constant and 



6(8) - 

For such u>(8) using (jnj) one has 



1 for 8 > 0, 
for 8 < 0, 



a(t) = ^ r 7 " -o e~li £t d8, (15) 



27T Jo rg _ £0)2 _|_ (2±\2 



2 



where 

i i r 7° 



N 2tt J ( g _ g oy + (7±v 



— d8. (16) 



Formula (fT5l) leads to the result 



. 7 o 
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I 2nl "'V h " * ' 2'* >", 
+ (-l) ft (-'(£S-^) t )]}, 



where E\(x) denotes the integral-exponential function [T3| [T4]. 
Using ( |T5i) or ( jTTl) one easily finds that 

ift^ = (£j-i 7 J)a(t) + Aa(*), 



where 



7T 2 



^(-i(«S + |7»«). 



So, 



where 



and 



hS) = ih^^^hl + Ah,(t), 



h° = f° - - a, 

'V — °<j> 2 ^ 



M4t) = ■ Aa[t) 



a(t) 

Using the asymptotic expansion of Ei(x) [14J, 
where | argz| < §7r, one finds 

— ih° + 

*(*)Uo * Ne * * 

+ ( j 2vr |*5|* U 
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and 



Aa(t)\ 



t^oo 



Nil (K¥ + 



2n \h$\* U. 

These two last asymptotic expansions enable one to find (see 

Aa(t) 



AM*)I 



t— »oo 



a(t) 



h 



hi- % - -2 



£1 /ft\ 2 



- +... 



Thus, there is 



M*)l 



h 



£1 ,h\* 



t— >oo 



— I — 



t \hl\l\tJ 



T + 



for the considered case (fill) of u(£). 
From (EZD it follows that 



def 



as t — > oo, where £?° = ^(i)! t_>oo, and 
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t— ► OOV 



7 



(as t — > oo). 



For different states |0 2 ) one has 



poo coo 
C 01 C <fe 



-2 



<70 

C 02 



^ - £l 2 * 0, 



ft\ 2 

7 



and 



whereas in general 7^ 7^ 7^. 
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{KM 



t-+oo/' 



The other relations are the following 



Coo 
C *i 


CO 

c <f>l 


1 h° 1 

1 n <f> 2 \ 


2 


Coo 

°<t>2 






2 




CO 

°<t>2 


i + l 


4>2 




c <f>i 




01 



^ % (32) 

^2 



01 02 

COO _ COO ,, |2 — ,,[) [J 

°01 °<fe 'Vl |fa 02 

01 02 



7^ 



c-0 ILO 12 _ CO \ h I 
CO I CO ' 

C 01 C </>2 



(33) 



and 



1 ? (34) 



C> (U , ( + \\ \ ' rsfl 



It seems interesting that relations ( J32l) - ( J34l) do not depend on time t. 

Note that the following conclusion can be drawn from (1301) : For suitably 
long times t there must be 



1*2 < K-^l- (35) 

These suitable times can be estimated using relation (|24p . From (1241) one 
obtains 

K*)UJ « + ^j^^ + ••• • (36) 

Relations (1261) — (1351) become important for times t > t as , where t as denotes 
the time t at which contributions to ^(t)^^ 2 from the first exponential 
component in (1361) and from the second component proportional to 4- are 
comparable. So t as can be be found by considering the following relation 

e n ^ 4^\h^¥- (37) 
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Assuming that the right hand side is equal to the left hand side in the above 
relation one gets a transcendental equation. Exact solutions of such an equa- 
tion can be expressed by means of the Lambert W function [15] . An asymp- 
totic solution of the equation obtained from the relation (1371) is relatively easy 
to find [16]. The very approximate asymptotic solution, t as , of this equation 



for (-$■) > 10 2 (in general for (-$) — > oo) has the form 
^ ~ 8,28 + 4 1n(%) 



h 



+ 2 In [8, 28 + 4 In (4)] + ••• ■ (38) 



3 Some generalizations 



To complete the analysis performed in the previous Section let us consider a 
more general case of a(t). Namely, let the asymptotic approximation to a(t) 
have the form 

"(*) ~ (39) 

t^oo k=0 

where A > and c& are complex numbers. Note that the asymptotic expan- 
sion for a(t) of this or a similar form one obtains for a wide class of densities 
of energy distribution u(£) [21 IH El El El E] ■ 
From the relation (139]) one concludes that 

9 -f ~ -X> + ^. (40) 

t^oo k=0 

Now let us take into account the relation (fTTl) . From this relation and 
relations (I3"9|) . it follows that 



M*) ~ T + F + f + - ■ (41) 

where c?i, c? 2 , c?3, . . . are complex numbers. This means that in the case of 
the asymptotic approximation to a(t) of the form ([3"9"]) the following property 
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holds, 

lim hJt) = 0. (42) 

t— >oo 

It seems to be important that results (|4"T|) and (1421) coincide with the 
results ( I27I) — (!30|) obtained for the density given by the formula (fT4|) . 
This means that general conclusion obtained for the other uj(S) defining 
unstable states should be similar to those following from (|2"T|) — (13"Uj) . 

4 Final remarks. 

Recently in many papers it was reported that some results of a detailed anal- 
ysis of the spectra of very distant astrophysical object could be explained 
by the variation of time of the fine structure constant a (191 HTJ [22], 123] - 
The analysis of the methods allowing for such an interpretation of the ob- 
served spectra shows that the conclusion that a varies in time was drawn 
by comparing the emission lines of two states separated by the fine-structure 
interactions and obtained by observing an astrophysical object with the cor- 
responding lines coming form the laboratory source [22] . In the light of the 
results obtained in Sec. 2 and 3 the question arises if such an interpretation 
of these astronomical observations is correct: For the observer the effect de- 
scribed by relation (l3"3l) and caused by long time properties of time evolution 
of unstable states, and the effect caused by a hypothetical variation in time 
of a, seems to be indistinguishable. Simply, there is a possibility that the 
interpretation of the mentioned results of astronomical results may not have 
taken into account all the available information. In the following we will try 
to shed some light on the above statement. 

So, let us consider a simple toy model. If one identifies energies £® with 
energy emitted by an electron jumping from the energy level rij to the energy 
level nk of the hydrogen atom in the nonrelativistic quantum mechanic then 
one finds 

= (4 -4), (n k <n j ), (43) 

n j k 

4 2 2 

where rij, n k = 1, 2, 3, . . . and R y = I! ^jS r = - ^ eC is the Rydberg constant, 
m e is the mass of the electron, c is the speed of light and a denotes the 
fine-structure constant. 

Now let these atoms interact with some external fields and let these fields 
excite the atoms forcing electrons to change their energy levels. Next, let 
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these excited hydrogen atoms emit the electromagnetic waves of the energies 
£° . = h v 1 ^ . (where v ® . denotes the frequency of the emitted wave) and let 
them move away from the observer with the velocity v. Then, the observer's 
measuring devices will show that the energies emitted by these moving atoms 
are £®' v . The energy £° emitted by this same source at rest is connected 

n jk 0,7 n jk J 

with £®' v D y the Doppler formula [18J 



^r=«^ fc < ( 44 ) 

where k = ]— — and 3 — -, and 

£°>* T £°>« = «(£° T £° ). (45) 
^From (jUJ) and it follows that 

CO,-!; CO ,.0, « ..0 

P o,« CO ' u ,o,« ,.o ' ^ v > 

c n 2k ^n 2k u ^2k n 2k 



where 82: v = hv® ,v , and similarly 

jh jk 



c0,v _ cO, v CO _ CO 

°rti fc °n 2 fc _ °ni fc °n 2 fc 

f0,» , o0,» ~~ CO _i_ fO 



(47) 



Now, let us assume for a moment that the to(£) given by the formula 
( TT4")) describes the considered sources with the sufficient accuracy. Next if 
one additionally assumes that in the distant past photons of energies, say 
£° , (j = 1,2), were emitted and that the sources of this emission were 
moving away from the observer with known constant velocity v, then at the 
present epoch this observer detects energies £^ = K>£™ k , where £™ k = f 
3? (h n . k (t)\. ■ If the moment of this emission was at suitably distant past 
then according to the results of Sec 2 (see (135]) ) it should appear that 

coo / coo '1^1 c0,v cO, v I (A%\ 

\°nik °n 2 k I ^ I °n lk °n 2k I 1 V* ) 

If this observer knows nothing about the effect described in Sec. 2 and 3 
and if he also finds that in addition 

coo cO 

-P + cT> (49) 

^n 2 k ^n 2 k 
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or 

Coo coo cO cO 

ri lk °n 2k i "ifc n 2k /rn\ 

Coo _i_ Coo i CO i CO ' ^ / 

then identifying ££° fc ' with (j = 1,2), and making reasonable assump- 
tions that - ( S^ k ' — £™' ) should also fulfil the relation of the type fH3|) . and 

the ratio = -0 s - or 0,0* should be the same as that at the present 

"2 fc n 2k n \k + n 2k 

epoch (see ). this observer, on the basis of property fl55|) . is forced to 
conclude that in the mentioned suitably distant past the Rydberg constant 
ftpast _ ^1^^ was smaller than at the present epoch. 

On the other hand, if the observer knows the effect described in the 
previous Sections then he could attribute the decrease of R y in the distant 
past to the long time behavior of the survival probability a(t) and resulted 
from this the decrease of the effective Hamiltonian h^t) for t —>■ oo (see (flOl) . 
( 1271) . ( 1281) and ( T4T1) ) rather than to an actual change in the fine structure 
constant. The conclusion that such an observation can be considered as a 
confirmation of the mentioned effect seems to be justified. Note that from 
( T46|) . fl33|) it follows that this effect seems to be distinguishable from the 
Doppler effect. 

Unfortunately, for the reasons mentioned in Sec. 1 it seems to be rather 
impossible to observe the discussed effect in laboratory tests. A similar con- 
clusion follows from the relations (|37j) . (|38j) . Nevertheless, this effect should 
be detectable by means of the observation of spectra of distant astrophysical 
object. 

As mentioned earlier some astrophysical observations suggest that rela- 
tions of type (T4gj) and (f5U|) take place. More precisely, these observations 
suggest that at very distant past the fine structure constant a was smaller 
than at the present epoch (for the recent data see [19]) which implies that 
R^ ast < R y . The known physical effects responsible for the observed red- 
shift are unable to explain such a result [20] • Therefore, the hypothesis that 
a varies in time, a = a(t), and that at the very distant past a{t) was smaller 
than a(0) (where a(0) denotes the value of a(t) at the present epoch) was 
formulated (see [19j [21], [22] and references one can find therein). The scale 
of changes of a is usually described by means of the quantity , where: 
Aa(t) = a(t) - a(0). 

One of the methods allowing to estimate consists in the analysis of 
the emission lines of two states separated by the fine-structure interactions. 
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In such to a very good accuracy (see formula (5) in |22j ) 



a 



L A 2 (*) + A a (t)J L A 2 (0)-A 1 (0)J' 



(51) 



where Xj(t), (j = 1,2), are wavelengths at time t. Using the relation £ n - k (t) 
^^y, (j = 1,2), one can rewrite this relation as follows 



Aa{t) 
a(0) 



+ 1 



x 



gn lfc (0)+gn 2fc (0) 
L£nu(0)-^ fc (0)J 



(52) 



Taking t — > oo in this formula and replacing £ n . k (t), (j = 1,2), by relations 



Aa(t) 



of type ffT2l . fl28|) . i. e., by £%° h , one obtains an expression for 
connecting the variation of a with the effect described in Sec. 2. Of course 



a(0) 



t—>oo 



the contribution of this effect into the possibly non-zero value of 



Aa(t) 
q(0) 



is only 



connected with the long time behavior of the survival amplitude a(t), that 
is with the long time properties of the real part of the effective Hamiltonian, 
£n h (t) = 3t(h njk (t)), described earlier and does not signify the change of a 
with time. 

Note that from (I3"3"l) and other results of Sec. 2 it follows that the right 
hand side of (j5"2|) should not be equal to 1 for t — > oo. So if one tries to 
estimate possible variations of a in time t using data obtained from astro- 
physical observation of a very distant cosmic object, one should take into 
account that a significant contribution to Aq obtained in this way can be a 
contribution generated by the long-time effect discussed in Sec 2. 

Unfortunately, the model considered in Sec. 2 does not reflect correctly all 
properties of the real physical system containing unstable states. Therefore, 
formula (1331) obtained within this model can not be considered as universally 
valid and it does not lead to the expression for ^ffl correctly describing real 
properties of very distant astronomical sources of electromagnetic radiation. 
The defect of this model is that some quantities calculated within it are 
divergent. Indeed from ([3]), (0 and ([T]) one finds that 

da(t) 



ih ■ 



dt 



H 



t=o 



(53) 



but inserting into (jSJ) the density to{£) given by (JHJ), relation fl53l yields 
(<p\H\(f)) = oo. For this same reason h^t) and Ah^(t) computed for u(£) 
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given by (TBI) are divergent at t = 0: /i«/,(= 0) = oo, Ah^it = 0) = oo. Nev- 
ertheless, taking into account that a very wide class of models of unstable 
states leads to a similar asymptotic long time behaviour of a(t) to this ob- 
tained for the model considered in Sec. 2, some general conclusions following 
from the results of Sec. 2 and 3 seem to deserve more attention. 

So, the following conclusion which can be drawn from (13"3"j) seems to hold: 
If one considers the same pairs of spectral lines (see formula (92) in [2T], or 
(5) in [22]) then the above discussed contribution into Aa calculated with the 
use of astrophysical data obtained for suitably distant two different sources, 
such that the second source is much older (i.e. much more distant) than 
the first one, can be the same for these different sources. What is more (see 
(155]) ). this contribution need not depend on time t but it should depend on 
such parameters describing the excited energy levels of atoms emitting the 
registered electromagnetic radiation and corresponding to measured spectral 
lines as energies 8®. k and decay rate ^° n . k (or lifetimes). This means that esti- 
mations of Aa performed by means of relations of type ([5"T1 - OK), (|52p and 
using the astrophysical data obtained for the given cosmic source emitting 
the electromagnetic radiation can lead to different values of Aa depending 
on the pairs of spectral lines used for calculations. If the hypothesis that a 
varies in time is true then it seems to be obvious that the above discussed 
contribution into Aa should be absent when one tries to estimate possible 
variations of a in time in laboratory test. Therefore one should expect the 
values of Aa obtained in laboratory tests [21J can differ from those obtained 
using astrophysical data. Simply the magnitude of possible variations of a 
obtained within the use of astrophysical data can be enhanced by the con- 
tribution of the effect discussed in this paper. 

All the above remarks concerning estimations of Aa using results of ob- 
servations of the spectra of distant astrophysical objects hold also when one 
tries to indicate a possible cosmological variation of the proton-electron mass 
ratio analyzing spectra of the same objects J2TJ [23]. Cosmic distances and 
other parameters computed from the observed redshift of very distant ob- 
jects emitting electromagnetic radiation |20j are calculated without taking 
into account the possible quantum long time energy redshift described in 
Sec. 2 and Sec. 3, so these distances as well as values of these parameters 
need not reflect real values of these parameters and thus need not reflect 
correctly real properties of the Universe. 

At this point it should probably be made clear that the discussion of this 
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Section and the results obtained in Sec. 2 and 3 are not meant to provide an 
alternative explanation of the source of the potential variation of a. Rather 
we point out there is a possibility to go beyond the standard interpretation 
of some results of the measurements of the scale of this potential variation. 

References 

[1] S. Krylov, V. A. Fock, Zh. Eksp. Teor. Fiz. 17, (1947), 93. 

[2] L. A. Khalfin, Zh. Eksp. Teor. Fiz. 33, (1957), 1371 [Sov. Phys. — JETP 
6, (1958), 1053]. 

[3] R. E. A. C. Paley, N. Wiener, Fourier transforms in the comlex domain, 
American Mathematical Society, New York, 1934. 

[4] M. L. Goldberger, K. M. Watson, Collision Theory, Willey, New York 
1964. 

[5] R. G. Newton, Scattering Theory of Waves and Particles, Springer, New 
York 1982. 

[6] L. Fonda, G. C. Ghirardii and A. Rimini, Rep. on Prog, in Phys. 41, 
(1978), 587. 

[7] A. Peres, Ann. Phys. 129, (1980), 33. 

[8] P. T. Greenland, Nature 335, (1988), 298. 

[9] D. G. Arbo, M. A. Castagnino, F. H. Gaioli and S. Iguri, Physica A 
227, (2000), 469. 

[10] K. Urbanowski, Phys. Rev. A 50, (1994), 2847. 

[11] J. M. Wessner, D. K. Andreson and R. T. Robiscoe, Phys. Rev. Lett. 
29, (1972), 1126. 

[12] E. B. Norman, S. B. Gazes, S. C. Crane and D. A. Bennet, Phys. Rev. 
Lett, bf 60, (1988), 2246. 

[13] K. M. Sluis, E. A. Gislason, Phys. Rev. A 43, (1991), 4581. 



15 



[14] Handbook of Mathematical Functions, Natl. Bur. Stand. Appl. Math. 
Ser. No 55, eds. M. Abramowitz nad I. A. Stegun (U.S. GPO, Washing- 
ton, D.C., 1964). 

[15] R. M. Corless, G. H. Gonct, D. E. G. Hare, D. J. Jeffrey and D. E. 
Khnut, Adv. Comput. Math. 5, (1996), 329. 

[16] F. W. J. Olver, Asymptotics and special functions, Academic Press, New 
York, 1974. 

[17] F. M. Dittes, H. L. Harney and A. Miiller, Phys. Rev. A 45, (1992), 
701. 

[18] W. Heitler, The Quantum Theory of Radiation, Oxford University Press, 
London 1954 and Dover Publications, New York 1984. 

[19] J. K. Webb, M. T. Murphy, V. V. Flaubaum, V. A. Dzuba, J. D. Barrow, 
C. W. Churchill, J. X. Prochaska and A. M. Wolfe, Phys. Rev. Lett. 87, 
(2001), 091301; J. Magueijo, J. D. Barrow, H. B. Sandvik, Phys. Lett. 
B 549, (2002), 284; M. T. Murphy, J. K. Webb, V. V. Flaubaum, 
Mon. Not. R. Astron. Soc. 345, (2003), 609; R. Srianand, H. Chand, P. 
Petijean and B. Brasil, Phys. Rev. Lett. 92, (2004), 121302; E. Peik, B. 
Lipphardt, H. Schneider and Chr. Tamm, Phys. Rev. Lett. 93 (2004), 
170801; C. J. A. P. Martins, A. Melchiorri, G. Rocha, R. Trotta, P.P. 
Avelino and P. T. P. Viana, Phys. Lett. B 585, (2004), 29; P. Tzanavaris, 
J. K. Webb, M. T. Murphy, V. V. Flaubaum and S. J. Curran, Phys. 
Rev. Lett. 95 (2005), 041301. 

[20] E. W. Kolb, M. S. Turner, The Early Universe, Addison-Wesley Publ. 
Co., 1993. 

[21] J.-P. Uzan, Rev. Mod. Phys. 75, (2003), 403. 
[22] J. N. Bahcall, Astrophys. J. 600, (2004), 520. 

[23] E.Reinhold, R.Buning, U. Hollenstein, A. Ivanchik, P. Petijean and W. 
Ubachs, Phys. Rev. Lett. 96, (2006), 151101. 



16 



